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Together with holonomy corrections, inverse volume terms should be taken into account when 
5-H , studying the primordial universe in loop quantum cosmology. We investigate how the tensor power 

spectrum is modified with respect to the standard general relativistic prediction by those semiclas- 
sical corrections. Depending on the values of the free parameters of the model, it is shown that the 
, spectrum can exhibit a very large deviation from its usual shape, in particular with a very red slope 

. and a strong running in the infrared limit. 

IT) , 

Q PACS numbers: 04.60.Pp, 04.60.Bc, 98.80.Cq, 98.80.C 
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INTRODUCTION 



A quantum theory of gravity is probably necessary to investigate situations where General Relativity (GR) breaks 
down. The early universe is a paradigmatic example of such a situation where the backward evolution of a classical 
. space-time inevitably comes to an end after a finite amount of time. Among the theories willing to reconcile the 
' Einstein gravity with quantum mechanics, Loop Quantum Gravity (LQG) is especially appealing as it is based on 
a nonperturbative quantization of 3-space geometry (see, e.g., [l[ and Q for an introduction). Loop Quantum 
Cosmology (LQC) is a finite, symmetry reduced model of LQG suitable for the study of the whole Universe as a 
simple physical system (see, e.g., Q). On the other hand, it is well known that the inflationary scenario is currently 
the favored model to describe the first stages of the evolution of the Universe (see, e.g., [J] for a recent review). 
Although still debated, it has received quite a lot of experimental confirmations, including from the WMAP 5- Years 
results Q, and solves most cosmological paradoxes. In this article, we consider the influence of LQC corrections 
to general relativity on the gravitational wave production during inflation. In this intricate framework, we assume 
the background to be described by the standard slow-roll inflationary scenario whereas LQC corrections are taken 
into account to compute the propagation of tensor modes. This approach is heuristically justified (to decouple the 
physical effects) and intrinsically plausible (as the LQC-driven superinflation can only be used to set the proper initial 
conditions to a standard inflationary stage if the horizon and flatness problems are both to be solved @). In 0,0, 
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holonomy corrections (due to the fact that loop quantization is based on exponentials of the connections, rather than 
direct connection components) were exhaustively considered. We now focus on the other fundamental LQC correction 
: the inverse volume -or density-operator- (due to terms in the Hamiltonian constraint which cannot be quantized 
directly but only after being reexpressed as a Poisson bracket not involving an inverse) . In the first section, the basic 
formalism is given together with the equations of motion derived in this framework. The second section deals with 
the definition of the power spectra and the question of initial conditions. Some analytical results are obtained in the 
third section. Finally, the fourth section explores the full parameter space with numerical investigations. 



I. FORMALISM AND EQUATIONS OF MOTION 



Loop Quantum Gravity relies on a Hamiltonian formulation of GR. The nonperturbative quantum effects associated 
with the LQG quantization procedure lead to effective semiclassical LQC equations. The associated hamiltonian 
constraint has been obtained by several different approaches 0, liol . [TH . The most important step in formulating 
LQG is to rewrite canonical gravity in terms of Ashtekar variables, which are the densitized triad Ef and the Ashtckar 
connection A\ where Ef = e"/|detef|, e"e| = q ab , q a b is the spatial metric, and A l a — T l a + K l a with T the spin 
connection and K the extrinsic curvature. The indices run from one to three. When written in Ashtekar variables, 
the matter Hamiltonian for a general space-time becomes [l3| 



det^l 2 v /|det£|, 

The terms involving inverse expressions cannot be straightforwardly quantized and must be regularized by a 
dedicated procedure [Til . [l5j . The expressions which result from this approach are rather complicated, and, in 
particular, are subject to a number of quantization ambiguity parameters. In principle, the spectrum for the inverse 
volume can be calculated exactly in isotropic LQC but the regularization leads to some ambiguities [lij]. The first 
important scale is set by ii = y^ylpi, where 7 is the Barbero-Immirzi parameter. It can be understood as the 
length above which space-time is roughly continuous and the inverse spectrum can be described by a continuous 
function. The second relevant scale is = ii y/ j/3, where j, which takes half integer values, is one of the ambiguity 
parameters. Above this scale the eigenvalues of the inverse operator follow the classical values, while they radically 
differ below I* . 

Loop Quantum Gravity introduces strong modifications to the dynamical equations in the semiclassical regime (i.e. 
when the scale factor a is such that ipi < a < £ r ). They come from the density operator [l7| 

with q = (a/£*) 2 . For the semiclassical universe, the quantum correction factor is given by [l8| 



D{q) = 5» r- 
1 2m 



f 



m + 2 



. -. ifn+2 I 1 im+2 



Q 



m + 1 



(Z +ir +i -sgn( g -i)i< ? -ir +i 



(2) 



m being an ambiguity parameter satisfying < m < 1. The cosmological dynamics with a field x is then governed 
by the following set of differential equations : 



H 2 _ 87TG 




D 



= 



(3) 



4HD 



V( X ) 



where H = a/a is the Hubble constant and a dot means differentiation according to the cosmic time. When a > £*, 
the Universe enters its classical regime and, in the limit a £*, leading to D(q) ~ 1, the usual Klein-Gordon equation 
is recovered for the inflaton field. In the semiclassical regime, a < it has, however, been shown that spectacular 
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modifications to the standard dynamics can be expected. For a scalar field driven dynamics, it seems that the field 
can naturally be excited up its self-interaction potential, setting the initial conditions for slow-roll inflation (see, e.g., 
@, HH, HI, HH, HH HE])- This is a very appealing feature of LQC as the requirement fa > "SMpi (for simple fa 2 
inflationary potential), imposed by observations, is rather difficult to set in the standard framework. Furthermore, a 
period of superinflation (H > 0) is expected to generically occur (see, e.g., [H, H3, Hl| ) , irrespectively of the detailed 
shape of the potential. 



Following the notation of [25|, [26 
the density operator is given by 



the equation of motion for tensor modes with quantum corrections coming from 



3H-~\h - 4v 2 ft = 0, 



(4) 



We note that such corrections are alternatively denoted by a in [29j . Because of quantum corrections encoded in the 
iS-term, the standard transformation from cosmic time t to conformal time rj (drj = dt/(a(t))) with the usual field 
redefinition <f> = ah, does not lead to a Schrodinger equation (see, e.g., Eq. (26) in 19]). In particular, there is an 
additional (anti)friction term given by —S/S. To recast the aforementionned differential equation into a Schrodinger- 
like equation, we switch from cosmic time to conformal time and redefine the field according to 

Finally, by decomposing the field 4> over its spatial Fourier modes, one obtains 



— + S 2 ( V )k 2 -V 1 ( V ) 



Mv) = o, 



(5) 



with a potential term given by 



„, , a" a' S' 3 fS'V IS" 



The prime should also be understood as a differentiation according to the conformal time. 
(S ~ 1), the potential term becomes V ss a" /a, which is the usual GR expression. 



(6) 



In the classical regime 



In addition to density-operator corrections, gravitational waves propagating in a FLRW background receive quantum 
corrections from holonomies [2i| . The influence of these LQC corrections has been studied in @, § ■ 

The values of ambiguity parameter n, which depends on the scheme adopted to quantize holonomies, range between 
— 1/2 and 0, though it was recently shown that n = —1/2 is favored [34[. Moreover, if n > —1/2, the holonomy 
corrections may become a major contribution to the effective mass of the gravitons at the end of inflation [|[ , leading 
to a rather intricate picture. In the case n = —1/2, the potential term reads: 



v 2 (v) 



a 
a 



2^7 



Ml, 



1-4 



H 2 (y) 
l 2 L 2 (v) 



(7) 



where L is the comoving size of a given patch. 

We will resist the temptation to combine V\ and V2 into a single equation which would describe both holonomy 
and inverse-volume corrections to the propagation of gravitational waves as our approach is to decouple, as much as 
possible, the different physical effects. This is, however, the next logical step in this study. Whatever the class of 
quantum correction considered, the equation of motion for primordial gravitons can be written as a Schrodingcr-likc 
equation: 



drj 2 



E k ( V )-V( V ) 



fk(v) = 0. 



(8) 



The difference — V(V)] can be seen as the effective squared-frequency 0^(77). The correspondence between the 

different terms is summarized in Table HI If E^(rj) > V{rf), the solution oscillates whereas it becomes a coherent sum 
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Density-operator 

LU1 1 L^l^ 11U1113 


Holonomy 


fk 
Mv) 
v(v) 


4>k = ahj\fS 


V>fc = a h 
k 2 

v 2 (v) 



TABLE I: Correspondence between the energylike and potential-like term. 



of an evanescent and an exponentially increasing wave if Ek(ri) < V(rj). Amplification of quantum fluctuations then 
arises when Ek(rf) < V(rf), i.e. when lo\ is negatively valued. 

From now on, we switch to our main working hypothesis: as in [f|, the background is supposed to be classical (i.e. 
described by the usual slow-roll inflationary picture) whereas the mode propagation is corrected by the inverse- volume 
LQC term. This makes sense as a phase of standard inflation is anyway mandatory after the supcrinflation regime. 
Furthermore, this allows us to understand in details the physical origin of the observed features. During slow-roll 
inflation, the scale factor is given by 

a(v) = 4 \v\~ 1 ~ e , 

with e the first parameter of the slow-roll expansion. The energy density is assumed to be related to the Hubble 
parameter via the standard Friedmann equation. The S operator is given by p9j 



S(q >1)«1 + Xq- K/2 , 



(9) 



with A and k two positive constants, not well constrained in homogeneous models (k = 2n with the notation of 
[29|). Setting ~ ip\ to remain consistent with our hypothesis of a full "standard" [i.e. classical) evolution of the 
background, one obtains the following energy and potential terms for density operator corrections: 



and, for holonomy corrections: 



Ek(v) 

v(v) 

E k (v) 



f 2A 
3e 



fc 2 , 

2^ 



3e 



A K (l + 2 e )(^i 



(_pi 

Co 



2v^7 3 (l + 4e) 



i«(l+e)-2 



(10) 

(11) 

(12) 
(13) 



II. ANALYTIC RESULTS 



The general equation of motion is far too complicated to be analytically solved in the general case. However, for 
some particular values of the parameters, the spectrum can be computed, at least in the infrared (IR) or ultraviolet 
(UV) limits. Those calculations are both useful by themselves and convenient to check the numerical results obtained 
in the next section. 

Throughout all this article, the convention of [33| is used for the normalization of initial states. The particle 
interpretation of the considered quantum field theory imposes the Wronskian W of the mode functions to be equal to 
—i. However, because we are working with rescaled quantities, the mode functions have to be normalized so that: 



16in 



The power spectrum then reads: 



(14) 



Ik 3 



(15) 
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As S ~ 1 at the end of inflation, the power spectrum can be safely approximated by 



2k 3 



(16) 



A. k = 1 and e = 0: UV limit 

When k = 1 and e = 0, the UV limit of the power spectrum can be analytically computed. The equation of motion 
can be rewritten as 



f0 
dvj 2 



(l + 2ZH)fc 2 



2 Z 
if \q\ 



0, 



(17) 



with Z = \(£pi/£ ) K = \(lp\/£ ) in this case. 

It is useful to define two regions: region I is such that the potential term can be neglected and region II is such that 
the time-dependent term in the energy can be neglected. Obviously, region I corresponds to r\ <C t] c where ry c is the 
potential crossing time defined by E(r) c ) = V(r) c ). Region II corresponds to —Z^ 1 < jj < 0. In region I. the mode 
functions are given by a linear combination (LC) of Airy functions: 



/>j(r?) ~ A r Ai(x) + S/Bi(s), 



with 



x= -fc 2/3 (2Z)- 2 / 3 (l + 2ZM). 
The two coefficients arc determined by the standard Wronskian condition 



(2Zk 2 ) 



2\l/3 



2ilm 



A\B I 



. 16tt 



leading to the natural choice 



A, 



B, 



2ny/2 



1/3 : 

kV2Z ) Mpi 



2nV2 



[kV2Z\ Mpi 

In region II, the mode functions are given by a LC of Coulomb wavefunctions: 



(18) 
(19) 

(20) 

(21) 
(22) 



b H {r)) ~ AnF^Z^k, k |r/|) + B//Gi(Z/2fc, fc |r?|). 



(23) 



In the UV regime (i.e. fc — * oo), the potential crossing time is roughly given by rj c oc —1/k and region I extends 
to values of r\ higher than —1/Z. As a consequence, the regions overlap and the matching can be done for conformal 
times such that 

fc" 1 < M < z- 1 . 

This region corresponds to the time when the density-operator is close to unity until horizon crossing. In this 
overlapping region, the Coulomb functions take the following approximate form: 



with 



bu ~ [An sin {-kr\ + ip u ) + B u cos (-h] + ip u )] 

7T 7 e Z 



(24) 



Vii = "J 



2A: 



z 

arctan I — I , 
4fc , 
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where 7 e is the Euler's constant and ipu ~ — ir/2 in the UV limit. On the other hand, still in the UV regime, the 
Airy functions can be approximated by sine and cosine functions: 



lim Ai(x) 



lim Bi(x) = — 



7T^/ \x\ 
1 



( 2 I |3/2 T 

cos I - \x\' - - 



I 3 N - 4 



7ry |a;| 



(25) 
(26) 



Taking the limit a; — ► oo (as /c — > oo and 77 is in the overlapping region) and performing a Taylor expansion in Z\r}\, 
this leads to 



1 (2Z 



1/6 

6j ~ -^= I ^- ) {[A 7 cos ((pi) - B/ sin (<£/)] cos (-fcr?) + [Aj sin ((pi) + Bi cos (<£/)] sin (-kr])} , (27) 



with (pi = k/ZZ — 7r/4. 

After matching the solutions, one easily obtains the coefficients in region II: 



An = 



D 



11 



2iV2n 
M P iVk 
iA 



ii' 



(28) 
(29) 



where the UV limit of ipu has been used. The power spectrum is derived by taking the asymptotic limit for small 
arguments of Eq. (|23|) . This can be performed by noticing that in the UV regime, the Coulomb wave functions are 
well approximated by Bessel functions. This leads to 



e Pi \ 2 /8r(3/2)\ 2 



7) {-^) 



(30) 



which coincides with the power spectrum in GR. 



B. k(1 + e) = 2: full spectrum and IR & UV limits 

In this particular case, it is possible to calculate the full power spectrum and to derive explicit IR and UV limits. 
Some algebra is, however, necessary. 

The equation to be solved is the following: 



drj 2 



fc2 _2 + 3e_ A/t(1 + 2e) ^|i| 



fc = 0. 



(31) 



With 



= k' - A«(l + 2e) ( ) , 



(32) 
(33) 



this can be written as 



drf 



K 2 + 2Zk 2 \n\ 2 - 



= 0. 



(34) 



To solve this equation, we will express it as a General Confluent Equation whose general form is given by 
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a w 



3d +v + w*._£ 



dw 
dz 



A i 



to = 0, 



A 2 



and whose solution is 



CiM (a, 6, h(z)) + C 2 U (a, 6, h(z)) 



M, U being Kummcr functions and Cj being integration constants. 
To write Eq. ([34]) as Eq. (|36|) . Ai must vanish. We define 

h(z) = (3z a 



(35) 



(36) 



(37) 
(38) 



where a and f3 are constants. The requirement Ai = reads as 

2A + ab-a + l = 0. 

A 2 must then be equaled to the (Ek(r)) — V(t))) term in Eq. ([56]) . With our choice for ft, and /, A 2 reads: 

A(A-l + 6a-a + l) ( 2n( b \ OL 2 p 2 z a 



- \a^\a--) - 



Taking a = 2, this can be written as 



. A(A + 2b-2) ( b\ 
A 2 = — 5 - 4/9 o - - - /3V 



which equals Ek{rj) — V(?7) if z = Identifying the other terms leads to 



= \V\ 
= 2, 



4/3 (a- - 



A(A + 2fe-2) 
2A + 26 - 1 



-2-3e, 
0. 



This can be easily solved to obtain 



13/ 12 
A = -2 ± 2V 1 + -9 e 



±i/cV2Z 



X 2 



1 3 / 12 
2 T lV 1 + T eT W2Z 



To determine the signs, we require the solution to converge to the usual GR solution when the LQC 
vanishing. This means: 

Urn <j) k = 4~z[D x J v {kz) + D 2 Y v {kz)\ , 
z — >0 



(39) 



(40) 



(41) 



(42) 
(43) 
(44) 

(45) 

(46) 
(47) 



(48) 

(49) 
(50) 
(51) 
terms are 

(52) 
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with 



When Z — > 0, one obtains 



3/12 



(53) 



o= T- 



4iV2Z 

which tends to ^fioo. In the limit \a\ — > oo, Kummcr functions can be rewritten with Besscl functions: 



hm [M(a, b, h(z))] = T{b) { — ) J TV {kz), 



kz 



A+l/2 



(54) 



(55) 



and b can be reexpressed as a function of cither A or v. When plugged into the general solution with Kummer 
functions this leads to 



lim fa = ^[E 1 J T „{kz) + E 2 Y TU (kz)} [1 + (z 2 )] 

Z — >0 



(56) 



where Ei are constants expressed with V functions and the C, coefficients. To ensure the equality with Eq. (|52|) . the 
(+) sign must be chosen for b and the (— ) sign for A (the sign of (3 is not relevant). The solution to Eq. is 
therefore 



exp ( —ik\/§r] 2 ' 



1 _ 3 
" 2 2 



I 2 T 4V 9 4 fcV2Z 



^if 2 3 , 

,1 + -,/: 



12 

~9~ 



e, ikV2Zif 



iK 2 



/ i 3 / io 
I 2 4 V 9 4W2Z' 



12 



1+ g V 1+ g-e,*fcV2Z^ 



After performing a first-order Taylor expansion in e, the parameters are given by 



.4 



-2-e, 



iK 2 



5 e 

« = 7 + o + /=• 

4 2 4fc\/2Z 

To explicitly derive the power spectrum, the solution is rewritten as 



(57) 
(58) 

(59) 
(60) 

(61) 



1 



(v) = —a exp -ik\ -x 2 dM(a, b, ikV2Zx 2 ) + C 2 U(a, b, ikV2Zx 2 ) 



(62) 



where x = \r]\. The integration constants are determined by the Wronskian condition. To compute the Wronskian, 
we will set h = k\/2Zx 2 and focus on the remote past. The solution is 



b k {h) = (k^y' 2 V A / 2 e - 8 * [C l M(a,b,ih) + C 2 U(a,b,ih)}, 
~ (kV2Z^ 



A/2 U\Bts[-a] ( f Ji \ ( h 

— j^Yfi — \ ^ CX P ( + Tm[a] In (/i) J + /i exp I — i— — ilm[a] In (h) 



with 



A = e^^d 



r(6) 
r(a) 



/' 



= 7rIm[a]/2 



C 2 + Cie* 



» r(fe) - 
r(a- &)_ 



(63) 
(64) 

(65) 
(66) 
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In the remote pas, the Wronskian can be written as 



-i kV2Z 



|A| 2 -|mI 2 



And the Wronskian condition 



leads to 



VhW h {<t>l<f> k } 



Setting A = 0, this means that 



|A| 2 -H 2 



d = 



8tt 

A+l/2 ■ 



M|j ( kV2Z 



C 2 



Mpi 



kV2Z 



-A/2-1/4 



-irIm[a]/2 



To obtain the power spectrum, the limit ?/ — > is taken in Eq (I62D : 



0) ~ C 2 H- A r( ^_, 1) (ikVtZrf 2 ^' 



Co 



T(a) 

r(l + e) 



r(| + e ) 



ifcV2Z 



. 5 , (»fc^)" 



0(77) 



res + 1 + » A - ) 



The spectrum is finally given by 



(M ^|r(3/2 + e)| 2 (2Z) 



-e/2-3/4 



jt3/2-e exp (_ rf) 
\ 4fc v 2.Z / 



( 



ril + i 



('A' 2 



2 ■ 



4kV2Z , 

This establishes the full tensor power spectrum. In the IR and UV limits, Im(a) — > ±oo, so 



r(5 + £ 

4 2 4fcV2^y 



iX 2 \ 



27rexp I — 



4W2Z 



A'* 



4fcV2Z 



3/2+e 



After performing the Taylor expansion in e and taking the IR and UV values of AT 2 , one obtains 



(67) 
(68) 

(69) 

(70) 
(71) 

(72) 



(73) 



(74) 



V, 



(IR) 



V, 



(UV) 



toy ^ 2lr(3/2 + e ) j 

^ P1 \ 2 /2 3 +T (3/2 + e) 



[2^(1 



-2c 



2 c k 3 exp 



Z 



tp\/2Z(1 + e) 
2fc 



l + (3 + 5e)-2 +0(k- 4 ) 



(75) 
(76) 



It can easily be seen from those equations that important modifications to the usual GR picture are expected in 
the IR regime whereas the UV behavior is equivalent to the GR one up to O (k~ 2 ). In the IR range, one can both 
notice a very red spectrum and a very strong running of the index. 



The tilt of the spectrum is given by 
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and the running by 



A IR ) 



AW) _ 



irs/2Z{\ +e) 



2k 



97 

_ 2e _ ( 3 + 5e ) — 



(77) 
(78) 



Aw) _ 



nV2Z(l + e) 
2k 

4Z 

(3 + 5e)-. 



(79) 
(80) 



III. NUMERICAL RESULTS 

To explore the full parameter space beyond the particular case k(1 + e) = 2, there is probably no other way to go 
than to perform a full numerical investigation of the problem. To this aim, a dedicated code was developed. Initial 
and final conditions, defined respectively in the remote past and at the end of inflation, are however required to 
perform such a computation. 



A. Boundary conditions 

To define the initial states, we first perform the following transformation: 

\t]\=e x , <t> k =u{x)e x ' 2 . 
The differential equation (jSJ with E and V given by Eq. (fl0|) and Eq. (fTTj) now reads: 



(81) 



a u 
dx 2 



1 + 2A(^J e ^ 1+ ^\e 2x k 2 -\K(l + 2e)(^ \ e< 1+e > - | - 3e 



u = 0. 



(82) 



In the remote past, the new variable x tends to infinity and the squared frequency in the above Schrodingcr-likc 
equation becomes dominated by the term proportional to e K -i 1J r ( -) x + 2x k 2 . With this approximation and using y = 
[k(1 + e)x + 2x]/2, it becomes a Bessel equation. For r\ — > — oo, the mode functions can therefore be written as a LC 
of Hankcl functions of order 0: 



4>k(v -> -oo) 



-v 



A k H {z v ) + B k Hl(z v ) 



with 



2V2A 



n(l + e) + 2\e 



* fc|77| [K(1+e)+2]/2 



(83) 



(84) 



The amplitudes of the mode functions are determined by requiring the appropriate Wronskian condition as defined 
in Eq. (TH]). Using 



W z (Hl(z),H (z)) = — , 

7TZ 



it can be shown that the Wronskian is given by 



W n {^A) = -- [k(1 + e) + 2] \A k \ z - \B k f 
2nV2 



The natural choice 



A k = 



Bi 



[«(l + e)+2]Mpi' 



(85) 

(86) 
(87) 
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Modes wave number 



FIG. 1: Primordial tensor power spectrum as a function of the wavenumber. The solid line is the analytical calculation whereas 
filled circles are the numerical results. 



will therefore be made to fulfill the Wronskian condition. 

At the very end of inflation, one can notice that whatever the LQC correction considered, it becomes subdominant 
when compared with the GR term. As a consequence, the equation of motion becomes in this regime 

1 d 2 4> k a 



4>k drj 2 a 

The solution is therefore given by a growing and a decaying mode. Clearly, to estimate the power spectrum, the 
growing mode is the most important one: 

MV -»• 0) ~ C k a(rj). (89) 
From this last expression, one simply rewrites the power spectrum as a function of Ck ■ 

op 

V T (k) = — \C k \ 2 . (90) 



B. Results 



For the numerical investigations, the solution is first taken as Hankel functions of order 0, given by Eq. (|83p . in the 
limit 77 — > — 00. It is then numerically propagated according to Eq. j5]) with a variable step fourth-order Runge-Kutta 
code. When 77 — ► 0, the numerical results arc fitted with Eq. (|89|) and the resulting value of Ck is used to compute 
the power spectrum. Figure [1] displays, for k(1 + e) = 2, both the numerical results and the analytical calculation in 
the IR limit obtained with Eq. (|75|) . The excellent agreement illustrates the reliability of the approach. Our code has 
also been tested with LQC corrections set to zero: the numerical results are also in agreement with the analytically 
known prediction of standard slow-roll inflation. 

A few comments can be made about the general structure of Eq. ([5]) , with the energy and potential terms given by 
Eq. (|10p and Eq. (jlip . First, it should be noticed that in the limit rj — > —00, the potential is always subdominant when 
compared with the energy. This makes meaningful the assumption of "asymptotically free" initial states. However, 
in the remote past, the LQC correction terms are dominant, both in the energy and in the potential. This is why 
the vacuum structure in the limit 77 — > — 00 is far from obvious and the intuitive interpretation of the result is very 
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FIG. 2: Primordial tensor power spectrum as a function of the wavenumber for different values of the slow-roll parameter e. 




Modes wave number k 



FIG. 3: Primordial tensor power spectrum as a function of 
panel) and n (right panel). 
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wavenumber for different values of the LQC parameters A (left 



difficult. Furthermore, if k(1 + e) > 2 the potential diverges when rj — > — oo. This unusual feature remains harmless 
as the energy term is, whatever the parameters, much higher than the potential term. When 77—5-0, i.e. at the end 
of inflation, both the energy and the potential are dominated by the GR terms. 

Figures [2j Fig. [3] and Fig. [4] display the primordial tensor power spectrum as a function of the wavenumber for 
different values of the physical parameters. Our fiducial model is defined by e = 0, A = 1, K = 2, Hq = 0.1. We 
remind that e controls the slow-roll whereas A and k are LQC parameters defined by Eq. ([9]) and Hq is such that 
Hq = £pi/£q with a{rf) = £q \rj\ e . The main result of this article, which appears in all the figures, can easily be 
noticed: the spectrum in strongly infrared divergent due to LQC corrections. On the other hand, as expected, the 
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FIG. 4: Primordial tensor power spectrum as a function of the wavenumber for different values of the parameter Ho- 



UV limit coincides with the GR prediction. 

Fig. [2] aims at underlining the general evolution tendency of the spectrum as a function of the slow-roll parameter 
and the last value (e=0.1) is deliberately chosen above the observationally allowed range Although it cannot be 
easily seen on the plot due to the scale, the spectrum exhibits, as it should, a — 2e tilt in the UV limit. On Fig. [3J 
one can notice that, as expected, the higher the value of the LQC parameters, the stronger the deviation from GR. 
Finally, the case of Fig. [¥]is slightly more intricate as H is fundamentally a background parameters which, however, 
couples to LQC corrections via the A term in the IR regime. 

On Fig. the critical wavenumber k cr u is displayed as a function of the parameters. We define k cr i t so that 
V r r(k cr it) = t 2,limk^ooP r T{k) (except for a nonvanishing e, in which case the criterion is V^^ G \k cr it) = 2V^ iG \kcrit)) , 
as an indicator of the transition wavenumber between the LQC-dominated (k < k cr it) and the GR-dominatcd (fc > 
k C rit) regimes. Interestingly, it can be concluded from those plots that this critical wavenumber is highly dependent 
upon any LQC parameter. To some extent, this feature is more observationally relevant than the amplitude of the 
effect which is anyway huge in the infrared limit. The higher the LQC correction, the higher the critical wavenumber, 
the smaller the physical scales submitted to LQC corrections, and the easier the observation. Although probably 
fortuitous, it is worth noticing that the transition scale is nearly proportional to the energy scale of inflation. The 
dependence upon the first slow-roll parameter is very weak in the allowed range, making the predictions reliable from 
the viewpoint of a test of LQC. 



CONCLUSION 



The influence of holonomy corrections during slow-roll inflation was derived in VA |8[. This article follows the 
same approach but considers inverse-volume terms (complementing the approach of [3a ] which was performed in the 
framework of superinflation) . Both analytical results (for a few particular cases) and numerical results (sampling the 
full parameter space) were obtained. The general behavior is a very substantial deviation from GR in the k — > limit. 
This deviation affects the amplitude of the power spectrum and, more importantly, both the tilt and the running of 
the index. There are several ways in which this work should be developed. First, it would be welcome to include 
simultaneously holonomy and inverse-volume corrections in the same differential equation for primordial gravitons. 
Then, and more importantly, corrections to the background should also be taken into account. Although several 
articles are devoted to this point, no fully consistent numerical study of LQC propagation and background corrections 
is yet available. Finally it would be interesting to investigate LQC corrections to the scalar spectrum [3r3 |. This last 
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FIG. 5: Value of the critical wavemimber as a function of the inflationary parameters e (top-left) and Ho (top-right) and as a 
function of the LQC parameters A (bottom-left) and k (bottom-right). 



point is probably the most promising one from the observational viewpoint. 
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